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PRIORITY INFORMATION 

The present application claims priority under 35 U.S.C. § 1 19(e) to U.S. 
Provisional application number 60/216,037 filed July 3, 2001, the entire contents of 
which are herein expressly incorporated by reference. 
5 TECHNICAL FIELD OF THE INVENTION 

The present invention is directed to path selection in networks. In particular, 
the present invention involves selecting a path that meets user defined criteria, 
through the use of a method that employs a Lagrange relaxation variable. 



BACKGROUND 

10 Modern telecommunication networks carry different types of traffic having 

different Quality of Service (QoS) requirements. QoS requirements can be defined 
using various parameters for a desired communication path. These parameters 
include bandwidth, delay, jitter and loss probability, among others. The ability to 
fulfill QoS requirements depends on the actual path selected in the network for 

15 carrying the traffic. Path selection therefore has to consider application dependent 
QoS requirements associated with the traffic. Without an efficient QoS routing 
algorithm to select the path that meets the specified QoS requirements, networks may 
fail to find a path and reject a request for a connection, even though there might be an 
appropriate path available to successfully establish the connection. In addition to the 

20 use of QoS requirements in path selection, the selection of a path with a minimum 

cost, which is also referred to as "administrative weight," is desirable, for example, to 
minimize the overall load on the network of the new path. 
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FIG. 1 is a view of a modern telecommunications environment. In FIG. 1, 
homes, businesses, schools, government institutions, hospitals, etc., all connect into 
the telecommunications network, seeking to transfer data, text, graphic images 
(including scanned images or computer generated graphics), audio (including voice, 
5 music, computer generated sounds and other natural sounds (animals, the ocean, a 
rockslide, etc.)), video, and combinations thereof. Coaxial cable, fiber-optic cable, 
fiber & coax, twisted copper lines, RF (including line-of-sight (LOS)) transmissions, 
and satellite transmission systems), and other means of connectivity link these 
institutions to each other. 

10 FIG. 2 illustrates various nodes, and links between these nodes in an 

exemplary communication system. The communication system includes nodes S, A, 
B, C, D, X, Y, Z, and T, and the various links between these nodes. Nodes, in 
practical applications, are telephone switching stations, routers, receivers and 
transmitters, or any electronic system capable of receiving a signal in, and 

1 5 transmitting a signal out. Links are the means of connecting the practical electrical 
nodes; LOS RF transmissions, fiber-optic cables, etc., as discussed above. 

A link represents the physical connection between two nodes, with no other 
node between the two nodes. A path is a set of individual links that moves a signal 
from its origin, S, to its destination, T (otherwise shown as "S-T"). For example, path 

20 p (S-A-Y-T), path q (S-B-Z-T) and path r (S-C-X-T) are composed of different sets of 
links, yet still represent the path S-T. Each link and node has specifications which 
describe the ability to pass data through the node or link, and what effect there might 
be upon the data. For example, there are delays (e.g., queuing and propagation 
delays), costs (e.g., hop count, monetary cost, or some measure of a link's capacity), 
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and data error rates associated with each node and link. Thus, the goal is to create a 
set of nodes and links that meets or exceeds the QoS requirements for that path. 

As an example of a typical communications network configuration, suppose, 
in FIG. 2, node S is an individual's home, and a call is being placed across the 
5 continental United States. The telephone call emanates from the residential home 
(node S) as twisted copper wire (link S-B), then is transmitted across the country 
(from node B) via a satellite transmission system (link B-Y). At the receiving end 
(node Y), the call is converted back from RF to twisted pair (link Y-T), which is the 
medium that enters the recipient's home or place of business (node T). It will be 
1 0 recognized and appreciated by those skilled in the art, that a separate line which 

emanates from S might be a CATV coax line (link S-A), or fiber optic line (link S-C). 
Depending on the application, the individual at node S could have a choice of which 
link to route a signal. Then, a QoS requirement might be established which is path 
dependent. 

15 It will be recognized that although the number of variations of links and nodes 

which can establish a path are potentially unlimited, in practice, the goal is to 
minimize signal degradation and maximize correct transmission of data. Thus, the 
need to optimize links between nodes and ultimately, the path S-T. 

The parameters that are used to satisfy QoS requirements can be either 

20 cumulative or non-cumulative. Non-cumulative parameters are those that do not add 
together or "accumulate" by the addition of subsequent nodes and links. However, it 
might be the case that the quality of the overall path is limited by the link or node 
with the worst non-cumulative parameters. Examples of non-cumulative parameters 
are bandwidth and access rights. The non-cumulative parameters can be easily 
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bandled by simply ignoring the links that do not fulfill the QoS requirements during 
the path selection process. Cumulative parameters can be either additive (e.g., delay, 
jitter and administrative weight) or multiplicative (e.g., loss probability). The 
simultaneous presence of more than one cumulative parameter is difficult to handle. 
5 The problem of finding a path which fulfills two simultaneous cumulative 

parameters or finding a minimum cost path which satisfies only one cumulative 
parameter is designated "NP-Hard." "NP" means "Non-deterministic Polynomial." 
See, for example, the book "Combinatorial Optimization" by W. J. Cook, W. H. 
Cunningham, W. R. Pulleyblank, and A. Schrijver (John Wiley & Sons, 1998) for 

10 more details, the contents of which are incorporated herein by reference. 

There are two classes of methods used to solve NP-Hard problems. Methods 
of the first class, referred to as "Optimal-Solution" methods, guarantee an optimal 
solution, but are complex and relatively slow. The second class is referred to as 
"Approximation" methods, and use a heuristic (trial and error) or approximation 

1 5 technique to solve the NP-Hard problem. Approximation methods usually run fast, 
but they do not result in the optimal solution. They do, however, provide a solution 
that is near-optimal. Approximation methods usually do not guarantee the quality of 
their results and hence they need to be thoroughly tested. Additionally, the solution 
for both classes of methods cannot be expressed as a polynomial. 

20 For both classes of methods, it is highly desirable to increase the theoretical 

guarantees on the running time of the method and the quality of the result of the 
method. With respect to Optimal Solution methods, one example is a "Breadth-First 
Search" method. This method determines the delay of all initial links, and then, 
building on that link, considers subsequent links, again, using only the one with the 
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minimal delay. For example, referring again to FIG. 2, to establish a path from node 
S to node T, the Breadth-First Search method would first determine the delay in links 
S-A, S-B and S-C. The method would then select the link with the least delay, e.g. S- 
C, and then determine which links out of C-B, C-X, C-Y and C-Z have the smallest 
5 delay. The delay for each link is determined and is added to the previous links, to 
maintain a running total. This continues until either the destination is reached, the 
delay constraint is reached, or there is no possibility to improve the path. However, 
the "Breadth-First Search" method is guaranteed to be fast, but can find a minimal 
path for only one cost (delay in this example). The Dijkstra algorithm is an example 

10 of a sophisticated Breadth-First search method. 

There are several well known methods of the Approximation class of methods. 
One example is "Fallback Routing" which tries to find an appropriate path which 
meets all the QoS requirements. It tries one path after another, until a suitable path 
has been found, or no path is found that satisfies all the QoS requirements. This 

15 method is very simple, fast and always gives an appropriate solution if it exists. 
There is no guarantee, though, of finding the optimal path and no guarantee of the 
quality of the path found. 

Both classes of well known methods, Approximation and Optimal-Solution, 
have significant drawbacks. Optimal-Solution methods provide an optimal solution, 

20 but are complex to use and relatively slow. Approximation methods run much 

quicker, but provide medium quality results and some types do provide a guarantee of 
the quality of the result. However, the better the guarantee of result, the slower the 
Approximation method runs. Additionally, it has been shown that for NP-Hard 
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problems, no algorithm exists that determines a solution, the running time of which, 
can be expressed in polynomial form. 

Thus, there has been demonstrated a genuine need for a method which is 
simple and easy to use, runs relatively quickly, and presents high quality optimal 
5 solutions for the problem of selecting a path in a telecommunications network that 
meets specified QoS requirements. The method of the present invention fulfills the 
long-felt need of choosing effective paths in modern telecommunication networks. 



SUMMARY OF THE INVENTION 

The invention is a method of determining a path in a communications system 
10 that minimizes a cost function and satisfies an additional constraint, these constraints 
referred to as QoS requirements, thereby providing a solution to the QoS routing 
problem. The method determines the path in a manner that is both efficient in time, 
and effective in quality of result. 

It is an object of the invention to provide a method that finds an optimal path 
15 in terms of minimal cost, the running time of which can be predicted. 

It is an object of the invention to provide a method that finds an optimal path 
in terms of minimal cost, the running time of which can be predicted and is shown 
to be a polynomial expression. 

It is an object of the invention to provide a method that finds an optimal path 
20 in terms of minimal cost, that provides guarantees for each solution. 
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It is another object of the invention to provide a method that determines an 
optimal path, out of many paths, which minimizes cost and satisfies a delay 
constraint. 

It is another object of the invention to provide a method that determines an 
5 optimal path, out of many paths, which minimizes cost and satisfies a delay 
constraint, in an efficient and economical manner. 

It is another object of the invention to provide a method that determines an 
optimal path, out of many paths, which minimizes cost and satisfies a delay 
constraint, the running time of which can be predicted. 
10 It is another object of the invention to provide a method that determines an 

optimal path, out of many paths, which minimizes cost and satisfies a delay 
constraint, the running time of which can be predicted and expressed in a 
polynomial format. 

It is another object of the invention to provide a method that determines an 
15 optimal path, out of many paths, which minimizes cost and satisfies a delay 
constraint, and provides guarantees for each solution. 

It is still another object of the invention to provide a method that determines 
an optimal path that minimizes cost and satisfies multiple constraints. 

It is still another object of the invention to provide a method that determines 
20 an optimal path that minimizes cost and satisfies multiple constraints, in an efficient 
and economical manner. 
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It is still another object of the invention to provide a method that determines 
an optimal path that minimizes cost and satisfies multiple constraints, the running 
time of which can be predicted. 

It is still another object of the invention to provide a method that determines 
5 an optimal path that minimizes cost and satisfies multiple constraints, the running 
time of which can be predicted and expressed in a polynomial format. 

It is still another object of the invention to provide a method that determines 
an optimal path that minimizes cost and satisfies multiple constraints, and provides 
guarantees for each solution. 
10 In accordance with one embodiment of the invention, the above identified 

and other objectives are achieved by determining a first path which has a first path 
cost and a first path delay, the first path cost being a minimum cost, determining a 
second path which has a second path cost and a second path delay, the second path 
being a minimum delay, calculating a Lagrange relaxation variable, determining a 
1 5 third path that minimizes a first modified cost function, calculating a second 

modified cost function for the first path, and determining the minimum cost and 
minimum delay path according to the first and second modified cost function. 



BRIEF DESCRIPTION OF THE FIGURES 

The novel features believed characteristic of the invention are set forth in the 
20 appended claims. The invention itself, however, as well as other features and 

advantages thereof will be best understood by reference to the detailed description 
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of the specific embodiments which follows, when read in conjunction with the 
accompanying drawings. 

FIG. 1 is an illustration of a modern telecommunications environment; 

FIG. 2 is an illustration of the switching path component of the 
5 telecommunications environment and demonstrates the nodes and set of links in a 
communication system; 

FIGS. 3 A and 3B illustrate a method for selecting a path in a 
communications system according to the preferred embodiment of the invention; 

FIGS. 4 A and 4B respectively illustrate a network and the Lagrange function 
10 of the network in which several paths have the same modified cost function at an 
optimal Lagrange relaxation variable; 

FIGS. 5 A and 5B respectively illustrate a network and the Lagrange function 
of the network in which an optimal path has a higher modified cost than a sub- 
optimal path at an optimal Lagrange relaxation variable; 
15 FIGS. 6A and 6B respectively illustrate the Lagrange function of a network 

using linear cost and using square cost; 

FIG. 7 illustrates a method for selecting a path in a communications system 
according to the preferred embodiment of the invention, with a first additional 
feature; 

20 FIG. 8 illustrates a method for selecting a path in a communication system, 

according to an alternate embodiment of the invention; 

FIG. 9 illustrates a graph which compares the average cost versus delay 
constraint for the present invention and other path finding algorithms; 
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FIG. 10 illustrates a graph which compares the average delay versus delay 
constraint for the present invention and other path-finding algorithms; 

FIG. 1 1 illustrates a graph which compares the routing steps versus delay 
constraint for the present invention and other path finding algorithms; 
5 FIG. 12 illustrates a graph which compares the routing steps versus delay 

constraint for various embodiments of the present invention and other path finding 
algorithms; and 

FIG. 13 illustrates a graph which compares the average cost versus delay 
constraint for various embodiments of the present invention and other path finding 
1 0 algorithms. 



DETAILED DESCRIPTION OF THE INVENTION 

The preferred embodiment of the invention solves the problem of determining 
an optimal path, in terms of one or more constraints, in a communications system, in a 
minimal amount of time. The present invention combines the well known Dijkstra 

15 algorithm (described below) with use of the Lagrange relaxation technique. 

Conceptually, any communication system can be thought as a collection of nodes and 
links. A collection of nodes is a "set" of nodes, and a collection of links is a "set" of 
links. A cost function can be written which describes cost as a function of the set of 
links. Recall that "cost" is simply any criteria that specifies a QoS for a link or node. 

20 There are constraint functions d u d 2 ... d k that describe the different constraints on 
costs that might be placed on a path in the network. The constraint functions d x , d 2 . . 
. d k have maximum constraints or limits D l5 D 2 ... D k . The task, therefore, is to find 
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a path/? from s to t (S-T) such that d,(p) < D i} (i.e. the constraint value of interest for 
path/? is less than the maximum allowable value for that constraint) for all i = 1 ... k 
and the cost of path p is minimal with respect to the cost function c among the paths 
satisfying the constraints. 
5 For purposes of illustration, path p, from S-T, might be the path consisting of 

the following links: S-A, A-X and X-T. Likewise, path q, from S-T, could consist of 
the following links: S-A, A-Y and Y-T. D is the maximum allowable delay, and d(p) 
is the delay associated with path p (assuming a path "x" exists), and c(p) is the cost 
associated with path p. The constraints d, describe the different constraints on paths 
1 0 (i.e., bandwidth, jitter, time-shift, etc.) that might be placed in the network. 

In one embodiment of the invention, there are only two constraints, cost and 
delay, to be minimized. The cost can be quantified in terms of any types of QoS 
parameters, including, but not limited to, number of hops, an additive monetary cost, 
or preference given by the administrator. The delay constraint is relaxed using the 
1 5 well-known Lagrange relaxation method. In the Lagrange relaxation method, the 
constraint is eliminated but it is built in the object function. That is, we eliminate 
looking directly to solve the problem of minimizing for the constraint in the paths, 
but, instead, indirectly solve for the constraint, by inserting the Lagrange variable 
(which represents the desired constraint) into the cost ("object") function, described 
20 above. In this way, what is originally very difficult to do directly, becomes much 
easier to perform indirectly. 

In performing the Lagrange relaxation method, a variable A>0 is introduced 
and for a fixed X, we are looking for an S-T path P x which approaches the following 
minimum: 
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L(A) = min{c(P) + A*d(P) - A,»D}; 
which can also be written as: 

L(A) = min{c A (P)} - A # D} . 

It is easy to find such path because X m D, i.e., the relaxed maximum delay 
constraint for any path, does not depend on the path P, so the task is to find a 
minimum cost path with respect to the modified cost function: 

c x (x) = c(x)+A.»d(x) 

where x is a path from S-T (and is different then P). 

It can be proved that L(l) is a lower bound of the original problem for any X > 
0. So, the best lower bound can be obtained by maximizing the function L(X). The 
term lower bound means the theoretical best solution. Let: 

L*=max{L(A): X > 0} 

and suppose X* maximizes the function L(X) [i.e. L(X*) = L*]. Then, it can 
also be shown that there is a path P* which minimizes the cost function c x and 
satisfies the delay constraint. Moreover, it can be seen that this path is the best that 
can be found by minimizing combined cost functions. 

FIGS. 3A and 3B illustrate a method for selecting a path in a communications 
system according to one embodiment of the invention. 



-13- 

The method presented below maximizes the function ~L(k) and finds the 
desired path P*. This method, the preferred embodiment of the invention, is very 
efficient in practice and its running time, i.e. number of steps, can be estimated in a 
polynomial expression. 
5 Initially, a path p between nodes S and T which has a minimal cost and a 

path q between nodes S and T which has a minimal delay are determined (step 301) 
through use of the Dijkstra algorithm. The Dijkstra algorithm is a very well known 
algorithm for finding a shortest path in a weighted directed graph. A directed graph 
is simply a set of nodes and links, as in FIG. 2. "Weighted" means that each link 

10 (and node) has a cost associated with it. Accordingly, the Dijkstra algorithm is 
performed once using cost of each link to find the minimal cost path and it is 
performed once using the delay of each link to find the minimal delay path. Next it 
is determined whether the delay of path p is less than or equal to the maximum 
delay constraint D (step 302). If the delay is less than or equal to the maximum 

15 delay constraint ("Yes" path out of decision step 302) then path p is selected since it 
is the minimal cost path which also satisfies the maximum delay constraint (step 
303). Then, the result is: 

k* = 0 and P* = p 

where X* is the A that maximizes the function L(A) and P* is the solution, 
20 i.e. the path that satisfies the delay constraints and minimizes the cost. 
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If, however, the delay of path p is greater than the maximum delay 
constraint ("No" path out of decision step 302), then it is determined whether the 
delay of path q, i.e., the minimal delay path, is greater than the maximum delay 
constraint D (step 304). If the delay of path q is greater than the maximum delay 
5 constraint D ("Yes" path out of decision step 304), then it is determined that there is 
no acceptable path between nodes s and t (step 305). 

If it is determined that the delay of path q is less than the maximum delay 
constraint D ("No" path out of decision step 304), then the Lagrange relaxation 
variable X, which is equal to the difference in cost between paths p and q, divided 
10 by the difference in delay between paths q and p, is calculated (step 306): 

^ = c(p) - c( g ) 
d{q) - d{p) 

Next a path r is selected which minimizes the modified cost function c x (r) 
(step 307). The path r is found through use of the Dijkstra algorithm. The 
modified cost function is defined as: 

15 c x (r) = c(r) + hd{r) 

where: 

c(r) is the cost of the r path; and 
d(r) is the delay of the r path. 
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The modified cost function c x (p) for path p is then calculated (step 308). It 
is then determined whether the modified cost function for path r, c x (r), is equal to 
the modified cost function for path p, cjp), (step 309). If the modified cost 
function for path r is equal to the modified cost function for path p ("Yes" path out 
5 of decision step 309), then it is determined that the value for the Lagrange 

relaxation variable X is a value which maximizes the function L(A) (the value of X 
determined in step 306) and that path q is a path which minimizes the modified cost 
function c^ (step 310). The method then has accomplished its goal of finding an 
appropriate path which is path q, as path p fails the delay requirement and path q 
10 fulfills it. 

If, however, it is determined that the modified cost function for path r is not 
equal to the cost function for path p ("No" path out of decision step 309), then it is 
determined whether the delay of path r, d(r), is greater than or equal to the 
maximum delay constraint D (step 311). If the delay of path r is less than the 

15 maximum delay constraint ("No" path out of decision step 311), then path q is 

replaced by path r (step 312), and the method reverts back to step 306. In step 306, 
a new X value is calculated, with the new path q (the path r just determined), and 
steps 307-309 again occur, with the new X value. 

If, however, the delay of path r is greater than or equal to the maximum 

20 delay constraint D ("Yes" path out of decision step 311), then path p is replaced by 
path r (step 313) and the method reverts back to step 306. In step 306, a new X 
value is calculated, with the new path p (the path r just determined), and steps 307- 
309 again occur, with the new X value. 
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The Dijkstra algorithm is used in step 301 (to determine two S-T paths:/?, 
minimal cost, and q, minimal delay) and in step 307 (to determine an S-T path r 
which minimizes the modified cost function c x (r)). Each time the Dijkstra algorithm 
is employed, it is for a single nonnegative cost or lengthcostcost function only, i.e. the 
Dijkstra algorithm can find a shortest path for one cost function only. 

Since there are only a finite number of different paths in the network the 
present invention will determine an optimal Lagrange relaxation variable A in a finite 
number of steps. The optimal Lagrange relaxation variable, A opt , is that smallest X for 
which the delay constraint of its associated "r" path is at most the maximum allowable 
delay constraint, D. The running time of the method of FIGS. 3 A and 3B can be 
shown to be expressed as a polynomial, because this method is an approximation of 
the NP-Hard QoS routing problem, and not an absolute solution. To find the running 
time of the algorithm, the notation 0(x) is used, which is defined as follows: 

g(x) = 0(f(x)) 

if there exists a constant c, c>0, such that: 

g(x) < c-fix), 
for all x. 

With this notation, the number of iterations (i.e. the number of executions of 
the main loop of the algorithm, steps 306 to 313 of FIGS. 3A and 3B), can be shown 
to be: 



0(mlog 3 m) 
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where m equals the number of links in the directed graph (or communication 
system). The only real time consuming parts of the main loop (and so in the entire 
algorithm) are the "calls" to use the Dijkstra algorithm (once in each iteration, plus 
twice during the initialization of "p" and "q"). Since the running time of the Dijkstra 
algorithm is expressed as : 

0(wlog m) 

the entire running time is expressed as: 
0((>wlog 3 »!+2)(»rlog m) = 0(m 2 -log 4 m). 

It has been recognized that there are two situations in which an optimal 
Lagrange relaxation variable X does not produce an optimal path. These situations are 
illustrated in FIGS. 4A, 4B, 5A and 5B. 

Referring to FIG. 4A, there is illustrated a multinode network, with six nodes: 
A, B, C, D, E and F, and eight links: A-B, B-F, A-C, C-F, A-D, D-F, A-E and A-F. 
Each link has a cost and a delay, shown in the parentheses, wherein the first number is 
the cost and the second the delay. There are four paths from A-F: A-B-F, A-C-F, A- 
D-F, and A-E-F. 

Referring now to FIG. 4B, there is illustrated a diagram in which all the values 
of L(A) for each path is plotted against X values equal to 0, 1, and 2. The "y" axis 
represents L(X), and the "x" axis represents X. Recall that the equation for L(X) is 
given by: 



L(A) = c,(r)-A-D. 
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In FIG. 4B, the values for L(A) have been plotted against the corresponding X 
values, and it is shown that all the paths, from A-F, have the same c x at the optimal X 
value. Table I illustrates the calculations performed to obtain L(A) for the various 
5 paths. At the optimal X (A opt ), the paths all have the same modified cost function, 
c x (r). In this instance, X opt = 1 and c x (r), for X = 1, is 7, for all paths. 
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1 5 Since all of the paths illustrated in FIG. 4A have the same modified cost 

function value, the present invention may not be able to select the optimal path. 
Although the situation illustrated in figure 4A appears to be problematic, in real 
networks this situation will not likely occur, and hence, is not worthy of any special 
consideration. 

20 Referring to FIG. 5A, there is illustrated another multinode network, with four 

nodes, A, B, C and D, and six links, A-D, A-B, A-C, B-D, B-C and C-D. There are 
five paths from A-D: A-D, A-B-D, A-C-D, A-B-C-D, and A-C-B-D. In FIG. 5B, 
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there is illustrated a diagram in which all the values of L(A) for each path are plotted, 
as was done in FIG. 4B. 

Table II illustrates the calculations to determine L(A), as was done in FIG. 4A. 
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However, FIGS. 5 A and 5B illustrate the situation in which a optimal 
Lagrange relaxation variable X does not produce an optimal path. Specifically, the 
network of figure 5 A will produce a modified cost function which is higher than a 
sub-optimal path at the optimal A. 

To address the situation illustrated in FIGS. 5 A and 5B, the influence of the 
cost in the modified cost function can be increased. FIG. 6A illustrates the 
Lagrange function of a network with linear cost. If, as illustrated in FIG. 6B, the 
cost were raised to the second power in c x , the effect of the cost becomes stronger. 
Therefore, the chance of the optimal path having both the smallest c x at the optimal 
A increases. Since the delay remains linear, the difference between the modified 
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cost function and the delay constraint, i.e., c x -D, remains a linear function in A for 
all paths. Further, it does not violate the convergence of the algorithm. Referring 
to FIGS. 6 A and 6B, it will be noticed that the modifying of the modified cost 
function c x does not influence the slope of the lines (i.e., c r D on A.), and instead it 
just shifts the lines farther off each other. 

FIG. 7 illustrates a method for selecting a path in a communications system 
according to the preferred embodiment of the invention, with a first additional 
feature. 

In this embodiment of the invention, the method is made more efficient when 
a path is found from a node to many other nodes. This is important, because, in 
implementation practical routing methods provide QoS work in routers and their 
task is not only to find an appropriate path to a certain other node, but they have to 
find paths to a bulk of destinations. The most resource consuming part of the 
method described with respect to FIG. 7 is finding minimal S-T paths with respect 
to the cost function c x (i.e. paths that are inexpensive, or, as inexpensive as 
possible). However, when the Dijkstra algorithm is "run" or employed in a weighted, 
directed graph (as shown in FIG. 2, for example), it generates path information for all 
paths originating at node S, to all possible destinations, and provides the minimum 
cost path for each source-destination pair. 

Therefore, in this alternate embodiment of the invention, whenever the 
Dijkstra algorithm is used with the modified cost function c A , the entire output of the 
algorithm is stored. This does not consume excessive amounts of memory because 
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the resultant paths form a tree. For improved efficiency, the costs and the delay of 
these paths are also stored. The data can be stored in either a binary search tree or 
in a balanced search tree in order to access it easily. 

When the Lagrange relaxation variables A were calculated previously, a path 
5 r was determined which minimized the modified cost function c x , and which is 
known as c A (r). Each path r, has a delay, d(r), and cost, c(r), associated with it. 
The first step in an improvement for selecting a path in a communications system 
(step 701) is to determine whether there exists a previously calculated path r having 
a delay greater than or equal to the maximum allowable delay D: 

10 d(r) > D 

If there is such a path, or paths, then in step 702 ("Yes" decision of step 
701), the r path associated with the largest Lagrange relaxation variable I is chosen, 
and is designated P^). The Lagrange relaxation variable A associated with that r 
path is known as X 1 . If there are no such paths, i.e. all the previously determined r 
15 paths have a delay d(r) less than the maximum allowable delay D ("No" decision of 
step 701), then in step 703 path p is selected as path P(A,). 

Similarly, in step 704, the delay of the r paths is checked to determine 
whether there exists a previously calculated path r having a delay greater than the 
maximum allowable delay D: 
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d(r) < D 



If there is such a path, or paths, then in step 705 ("Yes" decision of step 
704), the r path associated with the smallest Lagrange relaxation variable X is 
chosen, and is designated P^. The Lagrange relaxation variable X associated with 
5 that r path is known as X 2 . If there are no such paths, i.e. all the previously 
determined r paths have a delay greater than or equal to the maximum allowable 
delay D ("No" decision of step 704), then in step 706 path q is selected as P^. 

Following these steps, the method of FIGS. 3 A and 3B is followed, using 
path p equal to P(l t ) and path q equal to P(A 2 ). 
10 In an additional embodiment of the invention, an improved stop condition is 

implemented. An important feature of the method is that the path q always satisfies 
the delay constraint and it is replaced only when another path satisfies the same 
delay constraint and the cost of the new path is better than the previous one. 

As discussed above path p is a least cost path between nodes S and T, and 
15 path q is a path between nodes S and T which has minimal delay. If it determined 
that the cost of path q is within a certain percentage of the cost of path p, significant 
processing resources can be saved by simply selecting path q for carrying the 
traffic. This can be referred to as a "stop" condition, and can be performed between 
steps 304 and 306 of FIG. 3 A. This improves the running time, even if the error 
20 factor is very small, because only the last shortest path is eliminated from being 
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calculated in most cases. There is no significant loss of quality, since finding 
additional paths is done to ensure finding the optimal path. 

FIG. 8 illustrates a method for selecting a path in a communication system, 
according to an alternate embodiment of the invention. 

In this embodiment of the invention, the method described above can be 
extended to the case of multiple constraints. Additional constraints can be packet 
loss, phase shift, jitter, etc. In this case the method cannot guarantee that the 
resulted path fulfills all requirements but it attempts to reach a good compromise. It 
is immaterial the classification of the constraints, since a multiplicative constraint 
can be transformed to an additive constraint by using a logarithm approach (i.e., 
log(a-b) = log a + log b). 

In this general case a new variable X t , is introduced for each constraint dj. In 
this case, "d" does not just signify delay, it denotes a different type of constraint. 
For example, d x could be delay, d 2 could be packet loss, d 3 could be phase shift, etc. 
The modified cost function c k (x) (for some path x) is defined as: 

c k (x) := c(x) + EMW 

and L(A) is a lower bound of the optimal solution for all X, defined as: 



k 

LCX) = mm{c k (p) - £ Xp.} 
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where A : = (A l5 A 2 , . . . . A k ), and k is the number of constraints. Since 
the value L(A) is a lower bound of the optimal solution for all A, the goal is to find a 
vector A* which maximizes the function L(A) in the same manner as when there was 
one constraint. In addition we want to find a c k minimal S-T path. 
5 Initially, a number of paths are labeled p Q ,p 1 ,p 2 .... p k (which is one 

greater than the number of constraints) and stored (step 801). The path p 0 is a least 
cost (c-minimal) path, and p { are d x minimal paths for all i = 1 ... k. For example, 
if the four constraints are cost, delay, jitter and bandwidth, then path pj is the 
minimum delay path, p 2 is the minimum jitter path, and p 3 is the maximum 
10 bandwidth path. 

Next, in step 802, the solution X for the set linear equations is determined: 

c i(p<>) = c a(p,); V/ = 1 . . . k 

Once all the X values have been determined, a path r is calculated (step 803) 
that minimizes a multi-term equation for the modified cost function, c x (r), which is 
15 defined as follows: 

c 4 (r) = c(r) + ViW + X 2 d 2 (r) + X 3 d 3 (r) + A^r) 

Having calculated a path r that minimizes the modified cost function for a 
path r, the modified cost function is determined for the path p 0 (step 804). Then it 
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is determined whether q(p 0 ) = c k (f) (step 805). If the modified cost function for 
path p 0 is equal to the modified cost path for path r ("Yes" path out of decision step 
805), then the value of A has been determined which maximizes the cost function 
L(A). Then A* = A and P* = p t where path p t is chosen for which the maximal 
exceeding of the constraints is minimized (step 806). This is because any or all p t 
determined can exceed one or more constraints. 

If, however, c k (p 0 ) * c k (r) ("No" path out of decision step 805), then path p i 
is replaced with path r, (i.e. the equation from step 802 now becomes c A (p 0 ) = 
c A (r)), for all i (step 807). Each time p t is replaced with r (for each different i), a 
new set of A' s is calculated (step 807). The path p { is chosen which, when path r 
was substituted for it, maximizes the values of: 

i=l 

As an example, suppose there are three additional constraints, besides cost c: 
delay d x , jitter d 2 and phase shift d v 

We then need four paths (k+ 1 paths); p 0 , p,, p 2 and p 3 . The equations to 
find the solution for the modified cost function and A are as follows: 

(a) C M = c <Po) + K d i(Po) + <M 2 (Po) + MsCPo) 

(b) c A (Pi) = c(p x ) + KdiiPi) + X 2 d 2 (p t ) + A 3 d 3 ( Pl ) 

(c) %(P 2 ) = c(Pt> + Kdiip 2 ) + A 2 ^(p 2 ) + A 3 d 3 (p 2 ) 
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(d) c k (p 3 ) = c(p 3 ) + Vifo) + Kd 2 (Ps) + V 3 (P 3 ) 



The following are the equations to solve: 

1. c(p 0 ) + X^ipo) + X 2 d 2 (p 0 ) + A 3 d 3 (p 0 ) = c(p,) + Vi<Pi) + X 2 d 2 (p^ + 

l 3 d 3 (p 1 ) 

2. c(p Q ) + X^tpo) + X 2 d z (p 0 ) + l 3 d 3 (p 0 ) = c(p z ) + X^ipz) + X 2 d 2 (p 2 ) + 
l 3 d 3 (p 2 ) 

3. c(p 0 ) + X^dM + X 2 d 2 (p 0 ) + A 3 d 3 (p 0 ) = c(p 3 ) + A 1 fi? 1 (p 3 ) + X 2 d 2 (pJ + 
l 3 d 3 (p 3 ) 



Then, X l X 2 and X 3 are calculated. 

An improved stop condition can be implemented for the method of FIG. 8. 
It can be shown that c x (r) = L(X) < L* < c A (p 0 ) at each iteration step, so the method 
can be stopped when these bounds are close enough together, similar to what was 
shown before, in FIGS. 3 A and 3B. The improved stop condition, which can be 
inserted in the method of FIG. 8 after step 804, determines how far apart c A (r) and 
c i(Po) are > anc * saves significant processing resources. 

For example, it can be determined to stop the search for a path that satisfies 
the above described conditions when the two modified cost functions differ by less 
than 1 % . It is understood that the value of 1 % is for illustrative purposes only, and 
that other values can be used, dependent upon the circumstances. 
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Now that exemplary embodiments of the present invention have been 
described, a simulation of the present invention versus other path finding techniques 
will be described in connection with FIGs. 9-14 to highlight the advantageous 
characteristics of the present invention. The present invention is identified in 
5 figures 9-14 with the acronym LARAC which stands for Lagrange Relaxation based 
Aggregated Cost. The results illustrated in figures 9-14 are based upon a simulation 
environment of 100 random networks with 40 nodes and an average node degree of 
4 (node degree is the number of links connected to a node), which reflects real 
network values. The cost values on the links varies from 1 to 15 with a uniform 

10 distribution and the propagation delay on the links is selected to resemble delay 
characteristics of a nationwide network in, e.g., the United States. The first range 
of the propagation delay (1-5 ms) represents short local links, the second range (5-8 
ms) represents longer local links, while the third range (20-30 ms) represents 
continental links. The ratio of long local links was configured to be 20 percent, 

15 while the ratio of continental links was configured to be 5 percent. The results of 
the measurements illustrated in FIGs. 9-14 have been averaged. 

FIG. 9 illustrates the average cost for a variety of delay constraints. As 
illustrated in FIG. 9, all of the algorithms managed to find a path that satisfies the 
delay constraint. However, in minimizing the other additive metric (i.e., cost) the 

20 algorithms produced different results. The present invention found almost the same 
path as the constrained Bellman-Ford (CBF) routing algorithm, while a Delay 
Constrained Routing (DCUR) algorithm found paths with higher cost. The worst 
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cost was achieved by the Fallback algorithm. TABLE III below illustrates the 
maximum worst case cost max 6 A and the average worst case cost 5^ . 



TABLE III 



Algorithm 


max 5 A 


average 5^ 


Present Invention 


102.9% 


101.3% 


DCUR 


110.9% 


104.7% 


Fallback 


117.4% 


108.7%^ 


Lower Bound of Present Invention 


99.27% 


99.7% 



From FIG. 9 it can been seen that the first section of the curves (D < 15 ms), 
10 the found paths have very low cost, and after this delay point the cost of the paths 
increase with the delay bound. This can be explained by the fact that it is 
impossible to find a path for all source-destination pairs with a small maximum 
delay constraint D, and the found paths are short, hence, the costs of them are 
small. As the delay bound increases, the path-finding algorithms find longer paths 
15 and therefore the average cost of the paths increases. After the point where a path 
can be found for all source-destination pairs, the average cost of the paths will 
decrease with the delay bound since the algorithms are able to find paths with lower 
costs and higher delay bounds. 

FIG. 10 illustrates the average delay versus the delay constraint for the 
20 various path finding algorithms. As illustrated in FIG. 10, the curve of the present 
invention is very similar to the curves of the CBF, Fallback and DCUR algorithms. 
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However, the present invention has a smaller average delay than all but the 
Fallback algorithm when the delay constraint falls in the range of approximately 15 
ms to 35 ms. 

FIG. 1 1 illustrates the average number of steps of the various algorithms 
5 versus the delay constraint. As illustrated in FIG. 1 1, the present invention has an 
increasing running time for small delay constraints. This is due to the fact that as 
the present invention finds more paths to more destinations, the number of Dijkstra 
executions increases. As the present invention reaches the delay constraint that 
allows paths to all the nodes in the network, the number of steps starts to decrease 

10 since minimum cost paths computed first can more likely satisfy the looser delay 

bounds, thereby eliminating the execution of the second Dijkstra algorithm to obtain 
the minimum delay path. This is the reason why the average number of steps for 
the present invention approximates the number of steps of a single Dijkstra 
algorithm for high delay values. 

15 As also illustrated in FIG. 11, the Fallback algorithm has a decreasing 

characteristic as well; the difference is that it is linear in the first section and equals 
the running time of two Dijkstra algorithms. After this first section, the number of 
steps decreases to the number of steps of one Dijkstra algorithm for similar reasons 
to that described above in connection with the present invention. The DCUR 

20 algorithm was implemented in a centralized manner for the simulations. 

Accordingly, the running time of it would be significant, i.e., 2n times the running 
time of a Dijkstra algorithm. Thus, in FIG. 11 an estimation is plotted for a 
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distributed implementation. The estimation is the running time of two Dijkstra 
algorithms. The CBF algorithm has a totally different characteristic. Since the 
algorithm's exit condition depends only on the delay constraint, it increases with the 
delay constraint. The last horizontal section of the graph illustrates that the 
algorithm has been stopped because all the paths from the heap have already been 
used. 

FIGS. 12 and 13 illustrate the performance of the present invention, the 
present invention with the storing step, the present invention with a stop condition 
of 20%, 40%, 80% and 100%, and the Fallback algorithm. As illustrated in FIG. 
12, the gain of the routing steps under the critical delay bound of 15 ms is up to 
39%. FIGs. 12 and 13 illustrate the average routing step of the algorithms and the 
average cost of the computed paths. As can be seen from FIG. 12, the running time 
greatly decreases as we increase the value of the stop condition (illustrated in FIGs. 
12-14 as maximal difference (MD)). TABLE IV below summarizes the 
improvement of the algorithms' running time with the help of average iteration 
number. 

TABLE IV 
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rllgUl 1111111 


Average iieraiion 
Number 


max 8 A 


average 6 A 


Present Invention 
(unmodified) 


7.94 


102.9% 


101.3% 


MD 20% 


5.80 


104.7% 


101.5% 


MD 40% 


4.34 


105.8% 


102.4% 


MD 80% 


2.92 


108.4% 


104.6% 


MD 100% 


2.57 


110.2% 


105.4% 



As illustrated in TABLE IV, the unmodified present invention had an 
average iteration number of 7.94, while the present invention with a 40% stop 
10 condition the number of iterations decreased to 4.34. 

FIG. 13 illustrates the average cost of the path. When the stop condition has 
a 40% value, the average cost inefficiency increased to 102.4% . Similarly, when 
the stop condition has a value of 80% the average cost inefficiency equals 104.6%. 
The gain on running time compared to the loss on the optimality of path costs seems 
15 worthy. Moreover, the present invention provides a good way of controlling the 
trade-off between optimality of the path and the running time of the algorithm. 

The embodiments described above are merely given as examples and it 
should be understood that the invention is not limited thereto. It is of course 
possible to embody the invention in specific forms other than those described 
20 without departing from the spirit of the invention. Further modifications and 

improvements which retain the basic underlying principles disclosed and claimed 
herein, are within the spirit and scope of this invention. 



